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Abstract 

In a recent article we returned to the study of asymptotically flat 
solutions of the vacuum Einstein equations with a rather unconven- 
tional point of view. The essential observation in that work was that 
from a given asymptotically flat vacuum space-time with a given Bondi 
shear, one can find a class of asymptotically shear-free (but, in general, 
twisting) null geodesic congruences where the class was uniquely given 
up to the arbitrary choice of a complex analytic 'world-line ' in a four- 
dimensional complex space. By imitating certain terms in the Weyl 
tensor that are found in the algebraically special type II metrics, this 
complex world-line could be made unique and given - or assigned - the 
physical meaning as the complex center of mass. Equations of motion 
for this case were found. 

The purpose of the present work is to extend those results to 
asymptotically flat solutions of the Einstein-Maxwell equations. Once 
again, in this case, we get a class of asymptotically shear-free null 
geodesic congruences depending on a complex world-line in the same 



1 



four- dimensional complex space. However in this case there will be, 
in general, two distinct but uniquely chosen world-lines. One of which 
can be assigned as the complex center-of- charge while the other could 
be called the complex center of mass. Rather than investigating the 
situation where there are two distinct complex world- lines, we study 
instead the special degenerate case where the two world-lines coincide, 
i.e., where there is a single unique world-line. This mimics the case 
of algebraically special Einstein-Maxwell fields where the degenerate 
principle null vector of the Weyl tensor coincides with a Maxwell prin- 
ciple null vector. Again we obtain equations of motion for this world- 
line - but explicitly found here only in an approximation. Though 
there are ambiguities in assigning physical meaning to different terms 
it appears as if reliance on the Kerr and charged Kerr metrics and 
classical electromagnetic radiation theory helps considerably in this 
identification. In addition, the resulting equations of motion appear 
to have many of the properties of a particle with intrinsic spin and an 
intrinsic magnetic dipole moment. At first order there is even the clas- 
sical radiation-reaction term now obtained without any use of 
the Lorentz force law but obtained directly from the asymptotic fields 
themselves. One even sees the possible suppression of the classical 
runaway solutions due to the radiation reaction force. 



1 Introduction 

In a recent articlejl] we returned, with a rather unconventional point of view, 
to the study of asymptotically flat solutions of the vacuum Einstein equa- 
tions . The main development in that work was the realization that for any 
given asymptotically flat space-time with a given Bondi asymptotic shear, 
one can find a class of asymptotically shear-free (but, in general, twisting) 
null geodesic congruences where the class was uniquely given up to an arbi- 
trary choice of a complex analytic 'world-line ' in a four- complex dimensional 
space. Furthermore, by mimicking some terms that are found in the Weyl 
tensor of the algebraically special type II metrics, this complex world-line 
could be chosen uniquely and given - or assigned - the physical meaning of 
a complex center of mass; the real part being related to an asymptotically 
defined center of mass and the imaginary part related to asymptotically de- 
fined intrinsic spin. Equations of motion for the real and imaginary parts 
were found to strongly resemble the Mathisson-Papapetrou equations of mo- 
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tion for a spinning particle. 

The purpose of the present work is to extend those results to asymptot- 
ically flat solutions of the Einstein-Maxwell equations. Once again, in this 
case, we get a class of asymptotically shear-free null geodesic congruences de- 
pending on a complex world-line in the same four-dimensional complex space. 
However in this case there will be, in general, two distinct but uniquely given 
world-lines. One, which comes from properties of the Maxwell field, can be 
assigned as the complex center-of-charge while the other arising from prop- 
erties of the Weyl tensor, can be associated with the complex center of mass. 
Rather than study the situation where there are two distinct complex world- 
lines, we study instead the special degenerate case where the two world-lines 
coincide, i.e., were there is a single unique world-line. This mimics the case 
of algebraically special Einstein-Maxwell fields where the degenerate princi- 
ple null vector of the Weyl tensor coincides with a Maxwell principle null 
vector. Again we obtain equations of motion for this world-line - found here 
only in a second-order approximation. Though there is not a normal or stan- 
dard method for the assignment of physical meaning to different terms or 
quantities, nevertheless tentative assignments can be made by a variety of 
methods. One can use the Kerr and charged Kerr metrics to identify spin 
and magnetic moments. From the equations of motion that appear to have 
many of the properties of a particle with intrinsic spin and an intrinsic mag- 
netic dipole moment {the Mathisson-Papapetrou equations}, one can try to 
identify asymptotically defined centers of mass and charge. One even has 
at first order, the classical radiation-reaction term f now obtained with- 
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out any use of the Lorentz force law but simply from the asymptotic fields 
themselves. A new development occurs when we go to a 3rd order term in 
the equations of motion: we find a counter term that appears to suppress 
or dampen the exponential run-away solutions associated with the radiation 
reaction force. This arises solely from the Bondi mass loss equation. 

One should think of this work as developing a generalization of the prop- 
erties of the algebraically special space-times in the sense that the term that 
is required here to vanish, is automatically vanishing (among many other 
terms) for all the algebraically special metrics. It was, in fact, an under- 
standing of the algebraically special metrics and their associated shear-free 
null congruence that led us to this construction of the asymptotically shear- 
free congruences and the unique complex world-line. 

The charged Robinson-Trautman metrics and the charged Kerr metrics 
with their properties are explicit examples of the ideas and construction given 



3 



here. 

In Section II, which is a summary of our earlier review[I|, we discuss 
certain properties of null infinity 3 + that are needed here. In Section III 
and IV we describe the variables that come from the Maxwell and Weyl 
tensors and live on 3 + and their evolution, i.e., the asymptotic Maxwell 
equations and Bianchi Identities. In Section V we describe how to find the 
family of asymptotically shear free congruences, i.e., how to determine a 
field of complex stereographic angles, and then how to narrow that 

family down to two specific congruences: one comes from the Maxwell tensor 
while the other comes from the Weyl tensor. In the present work we shall 
assume that we are dealing with the special case when these two world-line 
coincide. We thus find, in this case, that there is a unique complex world-line 
in H-space that determines the unique shear-free null geodesic congruence. 
Finally, in Section VI, we attempt to give physical meaning and significance 
to the geometric structures that have arisen in this work. 

In particular, as a great surprise to us, we find that the equations of mo- 
tion for this complex world-line are extremely close to the standard equations 
of motion for a charged particle with intrinsic spin. The real part appears 
to represent a description of the motion of the center of mass of the interior 
gravitating system and the imaginary part the evolution of the spin vector 
of the same system. 

There appears to us a genuine mystery here: Why should the motion 
in a complex parameter space be virtually the same as ordinary motion in 
space-time? 

It is unfortunate, but apparently unavoidable, that there is considerable 
duplication of the material presented here with that of our earlier work[lJ. 
Not to have had this duplication would, we believe, make the present work 
unintelligible. 

2 Review of 3 + : the Future Null Boundary of 
Space-Time 

2.1 3 + coordinates 

Future null infinity, often referred to as 3 + , is roughly speaking the set of 
endpoints of all future directed null geodesies 0. It is usually given the 
structure of S 2 xR, a line bundle over the sphere. This leads naturally to 
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the global stereographic coordinates, (CO) f° r the S 2 and with u labeling 
the cross sections or global slices of the bundle. From the point of view of 
the space-time, 3 + is a null surface, with the null generators given by ((, () 
= constant. 

There is a canonical slicing of 3 + , ub = constant, referred to as the Bondi 
slicing, that it is defined from an asymptotic symmetry inherited from the 
interior, called the Bondi-Metzner-Sachs group[3j. Any one Bondi slicing is 
related to any other by the supertranslation freedom; 

UB = u B + a((X). (1) 

with a(C,C) an arbitrary regular function on S 2 . 

All other arbitrary slicings, (u — r), often called NU slicings, are given 

by 

u b = X(t,(,0- (2) 
2.2 Null Tetrads on 3+ 

In addition to the choice of coordinate systems there is the freedom to chose 
the asymptotic null tetrad system, (l a ,m a ,m a ,n a ). Since 3 + is a fixed null 
surface with its own generators, with null tangent vectors, [say n°], it is 
natural to keep it fixed. The remaining tetrad freedom lies, essentially, in 
the choice of the other null leg, l a . Most often l a is chosen to be orthogonal 
to the Bondi, ub = constant, slices. In this case we will refer to a Bondi 
coordinate/tetrad system. Relative to such a system any other tetrad set, 
(l* a ,m* a ,n* a ) 1 is given by a null rotationjlj about n a , i.e., 

j*a = la + hm a + l m a + }) l n a^ ^ 

m* a = m a + bn a , 
n* a = n a , 

b = -L/r + 0(r- 2 ). 

The complex function L(ub, C->Oi the (complex stereographic) angle be- 
tween the null vectors, l* a and l a at each point on 3 + , is at this moment 
a completely arbitrary angle field but later will be dynamically determined. 
Depending on how L(ub, (, C) is chosen l* a might or might not be surface 
forming. We will abuse the language/notation and refer to the tetrad sys- 
tems associated with l* a as the twisting-type tetrads as distinct from a Bondi 
tetrad even when l* a is surface forming. 
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The function L(ub, C> C) an d its choice will later play the pivotal role in 
this work. It will be seen that the vector field l* a can be constructed by an 
appropriately chosen L(ub, C C) so that it is asymptotically shear-free. 



3 Quantities Defined on 3 + for Any Given In- 
terior Einstein-Maxwell Space-Time 

Using the spin-coefficient notation for the three complex components of the 
asymptotic Maxwell field and for the asymptotic Weyl tensor, we have, from 
the peeling theorem El, that: 

0o = f + 0(r~ 4 ) (4) 
0! = ^ + 0(r~ 3 ) (5) 
2 = ^ + 0(r- 2 ) (6) 



= ^ + 0(r- 6 ) (7) 

^ = ^f + 0(r- 5 ) (8) 

^ 2 = ^f + 0(r" 4 ) (9) 

^3 = ;f + 0(r- 3 ) (10) 

^4 = ^ + 0(r- 2 ) (11) 



where O , X , 2 , -0 O , ■0 1 , -0 2 , -03, and ip A are functions defined on 3 + . For a 
given space-time their explicit expressions depend on the choices of both 
coordinates and tetrads. Their evolution is determined by the asymptotic 
Einstein-Maxwell equations and the choice of characteristic data. 



We adopt the following notation: expression written in the Bondi coor- 
dinate/tetrad system will appear without a star,'*', e.g., ip^, while for Bondi 
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coordinates with a twisting tetrad a '*' will be used, e.g., ip^ . In the case of 
NU coordinates and twisted tetrad, a double star will be used, e.g. i/j^* ,etc. 

Using Bondi coordinates and tetrad, the free characteristic data is given 
only by the (complex) asymptotic shear, 

a = a(u B ,(,0, 

while in Bondi coordinates, but with a twisting-type tetrad the free functions 
are 

<7* = <7> B ,C,C) & L(U B ,C,C) 

with (j*(ub,(X) an d L(ub,(,C) carrying the same (redundant) information 
as did the Bondi o~(ub, C 0- 

In the case of NU coordinates and twisted tetrad, the free data is given 

by 

a** = a**(T,(X) & V(r,C,C) 

with 

V(r,CX) = du B /dr = d T G(r,C,C)- 

The important point is that in both cases, the Bondi coordinates with 
twisting tetrad and the NU coordinates and twisted tetrad, all the informa- 
tion that was in <j(ub,(X) is transferred to the new variables, (cx*& L) or 
(<x**& V). Later we will show that all the information can be shifted into an 
appropriately chosen L(ub,CiC)i {which is also a complex function on J + } 
with a vanishing a* and that in certain special cases (a pure 'electric' type 
a) all the information can be shifted into the V(r,(,C) with vanishing a**. 

The relationship between the (0°, 0? 0o), ("04, ^3, V'?) ^0) given in a 
Bondi tetrad and the (0°*, 0?*, 0°*), (ipf, ipf, ip?, ip* °) of a twisting 
tetrad is 

0° = 0S° + 2L0f + L 2 0; , (12) 

0? = 4>f + i4f, (13) 

02 = <f>*2°, (14) 

^ = + 4L< + 6L V 2 ° + 4L 3 ^ + LV:°, (15) 
^0 = + 3L^° + 3L 2 ^ + L 3 ^l°, (16) 
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1>% 



1>f. 



(17) 
(18) 
(19) 



These are used, in the next section, going from the Bondi version of the 
Bianchi identities to the twisting version. 

4 The Asymptotic Maxwell equations and the 
Bianchi Identities 

The Asymptotic Maxwell equations and the Bianchi Identities in a Bondi 
coordinate and tetrad system are: 

A. In a Bondi coordinate and tetrad system[5j: 
Maxwell Equations 



0U- + g0 u _ a0! 
< + 50° = 0. 





'2 = 







(20) 
(21) 



Bianchi Identities 



1>l 




(22) 
(23) 
(24) 
(25) 
(26) 
(27) 



—a 



5 a — 5 2 cr + era' — aa', 



where 



By introducing the definition of the mass aspect 




(28) 



$ = + a V + aa' 



(29) 
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we see that (1271) is equivalent to 

* = (30) 

In the following section we show how the function L(ub,(X) can be 
chosen so that the shear a*, of the twisting congruence vanishes. For the 
moment we just assume that a* = for the twisting tetrad. 

By using Eqs. (|T2| - ([T71|) we find that for a Bondi Coordinate/Twisting 
Tetrad system with a* = 0, the same set of equations, the Maxwell and 
Bianchi identities, can be rewritten as: 



B. In Bondi coordinates with a twisting tetrad and a* = : 
Maxwell Equations 



C + + 2L'0* U + L<PY = 0, (31) 



Bianchi Identities 



*°- 1 50* 2 + {L<f)* 2 y = 0. (32) 



r °- 




+ 3k(j)* ' 


-*o 

?2 


(33) 




= -a^ 2 °-L^°--3L-^ 


+ 2k(f)f. 


-*0 

?2 


(34) 






+ tyfi 


-*0 
>2 


(35) 




= 5<7 + La" 






(36) 




= —a" 






(37) 




= 3 a — 5 2 cr + oo — aa 










+2L~5a" - 2LW + ZV' 


- L 2 W, 




(38) 



where 

a = bL + LL" . (39) 
By using (JUJ in (J2HJ) we finds that 

$ = + 2LM + L 2 a ■ + 5 2 cf + aa . (40) 

so that (|38|) is again equivalent to 

^-1> = 0. (41) 
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The dot denotes derivative with respect to «b- 



In the very special case where the Bondi shear is of pure 'electric type', 
the 'twist' of the asymptotically shear- free null congruence l* a can be chosen 
to vanish and hence the l* a would be orthogonal to the cuts of some constant 
r slices of 3 + . The slices, given by r = T(u B , (, () or ub = X(t, £, £), yield 
the parametric relations, (see below, Eqs. (jT7jl . (|4*K1 ) and (|5T|) ) 

L(ub,(,0 = a (T) x(r,cc) 
V(u B ,(,0 = kx(t,c_,0=x' 

u B = X(r,C,C). 

To perform the coordinate transformation to (r, (, () {from (ub, C C)}j 
referred to as NU coordinates, of the Maxwell equations and Bianchi identi- 
ties, it is convenient to note that 

3?7 = 3 (r) 77 - —rf 

V = ^rf (42) 

where 5( T ) is the edth operator taking r constant, r] is a function with spin 
weight s and rj is the derivative of r\ with respect to r. Then we have: 

C. In NU coordinates and 'twisting' tetrad with a** = 0: 
Maxwell Equations 

V<Po°' + 5(r)[U 2 0r°] = 0, (43) 

0r o, +o (T )[u0r ] = o. (44) 

Bianchi Identities 



= -V<5 (T) ijT° - 4[ 5 (r) V^r° + 3kV0* o *% (45) 

^o/ = _ya (r) ^*° -3[5 (r) U]Vr + 2kV(f)l*%* 

^> = -Vb {T) r 3 *° -2[5 (T) y]C° + kV^*%*° 

^ = v-%ft T) V + 28 (T) V] - V- 2 [f {T) V][b (T) V] 
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We emphasize that we have not lost any generality when going to the 
Bondi/Twisting tetrad but when we change to the NU coordinates the equa- 
tions are valid only when the original Bondi shear ct(ub, C was pure 'elec- 
tric'. In other words all the information in an 'electric' shear can be put 
into the real V(ub,Ci 0- Though in the present work these equations in 
the NU coordinates will not be used, they are included here for future use 
with further applications, including an analysis of the (charged) Robinson- 
Trautman[2j equations. 

5 Dynamics on 3 + 

5.1 The asymptotic shear- free condition 

We first show that there are special choices of L(ub,(X) leading to the 
twisting tetrads for which the shear, er* = 0. To find this family of func- 
tions L(ub, C C); a differential equation must be solved where the freedom in 
the solution is four arbitrary complex analytic functions of a single complex 
variable, i.e., an arbitrary complex analytic curve in a four-complex dimen- 
sional parameter space. [Surprisingly, the parameter space is the well-studied 
H -space.] All the information in the original Bondi characteristic data, 
aiusXiOi W1 ll have been shifted to the L(ub,(X)- Given an L(ub,(X) 
with any one of these curves, we could go backwards to recover the original 
Bondi o~(ub, C 0- The dynamics lies in the unique determination of this curve 
from other considerations that are describe later. The construction of these 
quantities, the L(ub,(X), the complex curve, the I*, etc., though done in a 
particular Bondi coordinate system, are invariant under supertranslations. 

We begin by observing that if we start with a Bondi coordinate/tetrad 
system with a given shear o~(ub, C C); then [after an unpleasant calculation 
the shear of the new null vector l* a , after the null rotation, Eqs.fJHJ), I s related 
to the old one by 

o-*(u B , CO = CO - 3L - LL . 

The function L(ub,(X) is then chosen so that the new shear vanishes, 
and hence, it must satisfy 

SL + LL- = <t(u b ,(,0- (46) 
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Remark 1 The special case, dL + LL = 0, played an important early role 
in the development of twistor theory, leading immediately to the relationship 
in flat space between shear-free null geodesic congruences and twistor theory. 
It leads us to conjecture that Eq. \46j) might play a role in some asymptotic 
form of twistor theory. 



Though Eq. lHEl) is non-linear with an arbitrary right-side and appears 
quite formidable, considerable understanding of it can be found by the fol- 
lowing procedure: 

Assume the existence of a complex analytic function of (m.b,C>C) {where 
C is allowed to be freed up from the complex conjugate of ( and where ub 
can take complex values}, 

r = T(ub, C,C) (47) 
that is invertible in the sense that it can, in principle, be written as 

u b = X(t,(,(). (48) 

Remark 2 We emphasize the obvious: since Ub is allowed to be complex, 
then the function X(t, (, Q is also complex. 

Then writing 

5T 

L = ~ (49) 
and using the implicit derivatives of Eq. P8|l , 

1 = X, T T (50) 
= 3 (r) X + X'tST, 

we obtain 

L = 5 (r) X (51) 

Again we mention that prime means the r derivative and 5( T ) means 5 with r 
held constant. Finally, with this implicit view of X and T, Eq. (fi6|) becomes 

5jV)X = a(X,C,C)- (52) 
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Remark 3 This equation had been derived earlier and was referred to as 
the 'good cut equation'^ In general it has a four complex dimensional 
solution space that has been named H-space. In its original form it had 
its origin in the search for complex null surfaces that were asymptotically 
shear-free. Here we are looking for real null geodesic congruences that are, 
in general, not surface forming, i.e., they have twist, but are asymptotically 
shear- freefTfy. 

Theorem 1 Note that using the reparametrization 

r ^ r * = F(r)=T*( Ms ,C,C) (53) 

leaves Eq. and the entire construction invariant. This fact will be used 
later for certain simplifications. 

Remark 4 We mention that if the shear in Eq. ^53\) . is pure electric [i.e., a (ub,(,() =3 2 S(ub, C? C) 
with S(ub, C? C) a real function] then the associated H-space is flat and has a 
real four- dimensional subspace that can be identified with Minkowski space. 

As just mentioned, the solutions to Eq. (|52j) depend on four complex pa- 
rameters, say z a , and can be summarized by 

u B = X(z a ,(,(). (54) 
However, since we are interested in solutions of the form, Eq. (j4*K|) . i.e., 

u b = X(t,(,(), (55) 

all we must do is chose in, if-space, an arbitrary complex world-line, z a = 
£ a (r), and substitute it into Eq. (|5H) to obtain the form (fHBTl . Furthermore, 
it can be seen that the solution, (IHH1) . can be written {with a special choice 
of H-space coordinatesjas 

u B = X(r, C, C) = e(r)l a (C C) + X l > 2 {r, (, C) 
with Xi>2 containing spherical harmonics, I > 2 and 

aU ' U " 2 li ' i + C c''i + cc'i + CC 



(56) 
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which has only the / = and 1, spherical harmonics. This result follows from 
the fact that 5^ in the good-cut equation annihilates the I = and 1 terms. 

The £ a (r) does get fed back into the higher harmonics via the a(X, £, £). 

Our solution to Eq. (|l6|) can now be expressed implicitly by 

L(u B ,C,() = 5 (T) X = ^(r)m a (C,C) + o (T )X i > 2 (r,C,C) (57) 
u B = X = e(T)l a ((,()+X l > 2 (T,(,(). (58) 

Remark 5 We emphasize that the complex world-line, £°(t), is not in phys- 
ical space but is in the parameter space, H-space. At this point we are not 
suggesting that there is anything profound about this observation. It is sim- 
ply there and whatever meaning it might have is obscure. This observation 
applies both to the vacuum Einstein equations as well as to the Einstein- 
Maxwell equations. Shortly we will see how to make unique choices of the 
world-line; in general the Maxwell field will have its own complex world-line 
while the Weyl tensor will determine a different world-line. We will consider 
here only the case where these two world-lines coincide. 

Remark 6 To invert the latter equation, $58\) . to obtain r = T(ub, C> C)> ^ 
in general, virtually impossible. There however is a relatively easy method, 
by iteration, to get approximate inversions to any accuracy. {See Appendix 
A} 

Remark 7 We have tacitly assumed throughout that the original asymptoti- 
cally flat space-time metric is real analytic or can be arbitrarily well approx- 
imated by a real analytic metric. 

If we expand £ a (r) in a Taylor series and regroup the terms with real 
coefficients and those with imaginary coefficients separately, we can write 

C(r) = G(r)+^(r). (59) 

This decomposition becomes important later. 
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Remark 8 The asymptotic twist, £(wb,C,C)j ( a measure of how far the 
vector field l* a is from being surface forming} is defined by 

2iS = 5L + LL -~5L- LL . (60) 

If L is obtained from a shear that is pure 'electric', then the C, a (r) can be 
chosen so that the asymptotic twist vanishes. 

It is from the fact that we can chose L(ub, C> C) as a solution of Eq. (j46|) . 
leading to a* = a** = 0, that is the justification for the form of Eqs. (E?T|) - (ITr]l . 



5.2 Determining the complex curve. 

The argument and method for the unique determination of the complex curve 
(and eventually the definitions of spin and center of mass motion) are not 
along conventional lines of thought. To try to clarify our argument, a brief 
detour might be worthwhile. 

It is well-known that for a static charge distribution, the electric dipole 
moment is given by the total charge times the center of charge position 
given in the 'static' Lorentz frame. In the case of a dynamic charge it is 
more difficult since the center of charge depends on the choice of Lorentz 
frame. Nevertheless by having, in any one frame, the center of charge at 
the coordinate origin, the electric dipole moment will vanish. In the case 
of a single charged particle moving on an arbitrary real world line, (the 
Lienard-Wiechert Maxwell field) the electric dipole moment vanishes when 
the coordinate origin follows the particle's motion or equivalently, the dipole 
is given by the particles displacement from the coordinate origin times the 
charge. The asymptotically defined dipole moment also vanishes when it is 
calculated (or extracted) from the / = 1 part of the coefficient of the r~ 3 
term of O = F ab l a m b , where l a is a tangent vector of the light-cones that are 
attached to the particles world-line. Analogously, it was shown[TT] that the 
magnetic dipole moment of a charged particle can be 'viewed' as arising from 
a charge moving in complex Minkowski space along a complex world-line. It 
is given by the charge times the imaginary displacement. {We emphasize 
again that the complex world-lines are a bookkeeping device and no claim is 
made that particles are really moving in complex space. The real effect of 
the complex world-line picture is to create a twisting real null congruence.} 
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If the light-cone from the complex world-line is followed to 3 + , (equivalent 
to introducing an appropriate asymptotic twist for the null vector l a that is 
used in O = F ab l a m b ) then both the asymptotic electric and magnetic dipoles 
vanish. This argument is given in considerable detail in the paper [ITj fT2]. 

The idea was to generalize this construction to GR. In linearized GR 
and in the exact case of the Kerr and charged Kerr metrics, the mass-dipole 
moment and the spin appear in the real and imaginary parts of the I = 1 har- 
monic of the coefficient of the r~ 4 part of the Weyl component, ip® , in a Bondi 
coordinate/tetrad system. This observation was extended to all asymptot- 
ically flat vacuum solutions by going to an asymptotically twisting tetrad, 
with an appropriately chosen iJ-space complex curve [see below], such that 
the new Weyl component, , has a vanishing I = 1 harmonic in the coeffi- 
cient of its r~ 4 part. {Note that for algebraically special type II space-times 
the entire quantity ip*® = 0, so by requiring that the coefficient of the lowest 
harmonic ,1 — 1, should vanish, we are imitating the type II properties as 
closely as possible.} We refer to this curve as the intrinsic complex center 
of mass world-line and, roughly speaking, it will identify the center of mass 
and the spin angular momentum from its real and imaginary parts times the 
mass. Effectively, by choosing the world-line so that the / = 1 harmonic van- 
ishes, we have shifted - in some sense - the 'origin' to the complex world-line. 

To apply this idea to our discussion of twisting null congruences for the 
Einstein-Maxwell case, we first treat the Maxwell equations. From the pair, 
(HH) and (JH3), that is 

0*o- + 9 0*° + 2L(j)f + L(f>l°- = 0, (61) 
0f + 50* 2 O + (L^ )- = 0, (62) 

we require that in (RUT) , the / = 1 part of 0q° should vanish. This deter- 
mines, in principle, the complex world-line, £ k (r), and leads to the definition 
of the complex center of charge world-line. {In practice we must use approx- 
imations for this determination.} In this work, as was mentioned earlier, we 
are making, though it is not necessary, the special assumption that the two 
complex world-lines, the complex center of charge world-line and the complex 
center of mass world-line [determined later] coincide. 

Remark 9 In what follows we use the recently introduced Tensor Spin-s 
Harmonics} 131/. Y"if] ., defined from the tensor products of the three basic Eu- 
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clidean vectors (c^m^mj). See Appendix C for some notation and gr-qc/0508028 
for details. In particular (Y^ = — q = — 2l i} = rrii, = mi). 

By using the freedom in the reparametrization of r given by (|5Tfl) and 
Eq. (f56|) in Eqs. ([5?]) and (f58|) and expanding, we obtain that (see appendix 
A for details) 

r « u B -C (ubMCQ+X^ubXX) (63) 
t(r) = (^V),^(r)) = ( v / 2r,^(r)) (64) 

e(r) « ^(«b) + irW^W^ (65) 
and hence from Eq. ljHTjl and 

z w 1 

we find 



= -^y^ + -y^. (66) 



£(«b, C, = lt k (u B ) + -{u B )e{u B )]Yl k + y> 2 . (67) 

Since 4>*\ and 02° have spin weight and —1 respectively, then they can 
be written in the following form 

0*° = c + C l Y^ + ... 

= GT^ 1 + ... (68) 

For our approximations, (essentially, perturbations off the Reissner-Nordstrom 
metric), we choose C as a zero order quantity, and C( u b), C % (u b ) and G 1 {ub) 
as first-order. We are considering only I = and 1 terms, i.e., the monopole 
and dipole. 

From Eqs. (JHTj) and (|68|) we find that 

50*° = -2C l Y^ + ... 

Ltf = {C-\e + ^JkCV] + ■j=e ijk eC>}Y* e + ... 
L4>f = {^>^^ j ] + ^a^ J K + ---> (69) 
where we have used the Eq.fJUHJ). 
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Therefore, working with the leading terms, we find that the condition 
[0o ]«=i = 0, is equivalent to: 

+C\t + ^mtC] + ^mC ci - (70) 

To see the meaning of this condition we go back to Eq. (JH21 )■ By using 
Eqs. (|H7j) and (fHHl) . a direct computation shows that: 

L$ = l^ j G j --^e ijk e^-Y 1 l} + ..., (71) 
where we have used the product decomposition 

Therefore, considering only the leading terms and the I = and I = 1 
harmonics we find, by using Eqs. (fTT| . that Eq. (|62|) requires 

C = (72) 
C f - = -G' + ^eyftfC?]-. (73) 

If we assume, for the moment, that C( T ) an d therefore C( u b) is known, then 
Eqs. (JZOJ) f72|) and (f73|) are a system of three differential equations for C, 
C % and G\ Allowing only second order terms of C( u )i the solution to this 
system is given by 

c = q ■ | ee- • .... 
c f = • .... 

G f = -q[t f + ^e ijf eer + .., (74) 
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where q is a constant. Therefore, the condition that 4>* should have no I = 1 
part, implies that the solution of the Maxwell equations up to I = 1 is given 

by 

0* 2 ° = -Q^+^jfeen-f 1 - (75) 

Observe that 02° can be written as 

<f>* 2 ° = c(P 2 = -D*~Y-\ (76) 

where 

Df(u B ) = q[i f {u B ) + -L=e ijf e(u B )e{u B )), (77) 

which, up to our order of approximation, is the complex dipole moment. This 
means that if we know the complex dipole moment then the complex cen- 
ter of charge can be obtained. This result had been obtained earlier|l2] for 
Maxwell fields in Minkowski space using the following procedure: first the 
asymptotic Maxwell equations in a Bondi coordinates/tetrad system were 
integrated assuming that the complex dipole moment is given and then by 
using the twisting tetrad version of the Maxwell field, (p* °, the C (r) was de- 
termined by requiring that the I = 1 part of (f)*® should vanish. 

Note that in above argument, we have not obtained any dynamics for the 
complex curve C( T )', a U we did was relate it to the electric and magnetic 
dipole moments of the Maxwell field. The dynamics is found by going to the 
Weyl tensor. 

Explicitly we study the asymptotic Bianchi identities under the condition 
that the I = 1 part of ipi° should vanish. We retain the assumption that the 
I = 1 of 0o° vanish with the same choice of C( T )i i- e -) that the two complex 
curves coincide. 

Starting with the set of equations (pjj| and (jHSJ), namely 

r o- = - L r 2 °- ~3L-r 2 ° + 2k<Pl% (78) 

r 2 °' = -5^ O -L^ o --2L-^ o + ^ o 0; (79) 
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we require, that in (f7H| . the I = 1 part of should vanish. This leads (via 
approximations) to differential equations for the £ a (r). Its evolution will be 
seen to be driven by the original Bondi shear and the Maxwell field. 

Before turning to the dynamics of £ a (r) we recall the definition of the 
Bondi mass-momentum and its relationship to the mass aspect ^ and its 
evolution. Eq. ([TH|l can be rewritten as 



= a a + k0* 2 % , (80) 
o = QL + LL' (81) 
= m - 2LM - L 2 W - 6 2 cf - aa. (82) 

<P*2° = =-vie + ^e ljf eeiY 1 - f 1 + (83) 

Using the definition[lJ[see Appendix B] that the Bondi four-momentum, 
P a , can be extracted from \I/ by 

P a (u B ) = (Mc, -P l ) = --^ / m a dS (84) 

07TG J 

then, by inversion, one can see that 

CO = - |FT» + ^> 2 (85) 

where ^> 2 contain only harmonics I > 2. Since \I/ is real it follows that both 
M and P % are real. 

The time evolution of P a (u) is found from Eq. (EDI) by integration over the 
sphere, i.e., 

P a {u B ) = J (aa + k<p* 2 %°)l a dS, (86) 

the Bondi mass/momentum loss. 

If ^*°, from Eq.dS2J, is now substituted into (f7H| we obtain after using 
Ea. pUl and regrouping, a truly ugly equation that can be written as 

^*o- + m + 3L y _ + k[L(f)* 2 % - 2<f>f 0°] = (NL), (87) 
with the non-linear terms, (NL), expressed by 

(NL) = Q[aa) + 2aSa + 3L5 2 a + L 3 a' + 3Lctct ' + 3L 2 5a" (88) 
+3L[2L5a + L 2 a ' + 5 2 a + ctct]. 



20 



The grouping has the non-linear higher order terms put into (NL) while 
the 'controllable' (lower order) terms are on the left side. It is from this 
equation and (|86|) we find the equations of motion. We require, as said 
earlier, that ?/>*° (and of course ip^ ' ) have no I = 1 harmonics. Since 5 3 a has 
spin s — 2, there is no I — 1 term and thus we see that Eq.pTjl reduces to 
the condition that 

[9* + V + k(L<f);% - 20*°0°)] i=1 = (NL) l=1 . (89) 

It is this equation that we must analyze, or more accurately approximate, by 
looking only at the leading spherical harmonic terms. From Eqs. f)85l ) and 
(JfiTj) . i.e., from 

*(«,C,C) 

L(u, (,() 
it follows that 



2V2G 



QG 



M - —I }i Yu ■ ■■■■ 



[e(u)+^e ijk e(u)e(u)}Y l \+... 



5^ = 12^%^ (90) 
GV 2G M[e + -^=e ijk tm k - ^meP^i + Y{1 > 2} 



t*0' 



7*0 



-far - ^AVfmi + rcvr»\ 

= -q 2 {f" - ^f^iXr r + -^e ijk fe}Yl k + Y{1 > 2}. (91) 

where we have used Eq. (|66| . Then, from Eq. (f89|) . keeping only the first and 
second order terms, we have from the I = 1 harmonic 



^f{f" - ^^''y - 2 te\}Yl k = 0, (92) 
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our basic equation which determines the Bondi momentum in terms of the 
complex world-line. 

Rescaling «b to the conventional Bondi coordinate 1 u c , by ub = ^ cu ci 
we have our basic equation, the Bondi momentum in terms of the complex 
world-line; 

-%{t'--^A^ry-2re]}- (93) 

footnote 1 : Many years ago, for reasons of symmetry and convenience the 
conventional radial coordinate r c was rescaled by a factor of a/2, i.e., 
r = y2r c , so that to keep a coordinate condition the retarded time 
coordinate u c was rescaled so that our ub = ^ u c- Here we will restore 
the conventional u c but also include the ordinary time units {i.e., with 
c^ 1} we take ub = ~^ cu c, so that derivatives d UB = -^d Uc . 

By using the decomposition 

and then iterating, i.e., by substituting into the right-side of Eq.fJHHl) the real 
linearized expression for P k , 



into (j93|) . we find that the Real Part of Eq. (|93|) is given by: 



-^,a^k j h - e, Cn + e,e,i + e R erh (94) 

while the Imaginary Part is 

M la 2 

= m# + Yc e ijk [&e R - ##] + m 

q 2 i j— i j— i.. 
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We can now obtain the equations of motion for £^ by taking the Ub 
derivative of Eq. (|M|) . i.e., 



3 . 1 



- d'G + r;e?r + ekf }•■ (96) 

and equate the right side to the Bondi mass, momentum loss, i.e., to Eq.lfHHj) 

P a (u B ) = J (a a + k(j>f^)l a dS. (97) 

k = (98) 

{Note that care must be exercised when using this equation since all the 
d UB derivatives must be replaced by the rescaled derivatives ^d UB .} 

In order to evaluate the integrand and the integral, we assume only the 
lowest gravitational harmonics, i.e., the quadrupole radiation and hence take 



a = o ij {u B )Y*+Y l > 3 , (99) 



where Y 2 \^ = fn^j- Then 



a-a = a^Y^Y^ + ...., (100) 



with Y 2kl = Y 2 \ y After a direct but lengthy computation P3| to evaluate the 
products of the harmonics, we find that 

a a = Irf '-n"- + + Y l > 2 . (101) 

Likewise we have that 

= yfr - ^e^re^l + Y> 2 . (102) 

Putting Eqs. (jl02jl and (llOljl into (f97l) the evolution for the mass and the 
three momentum is then given by 
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2a 2 —- r 

= -%[c R c R + era - ^ ,J ^ r - ( 103 ) 

• 2 • 2 

P fe - 1< 1 <: ~'"rJ- 10 - 7?U-~jl>- 

P = g^-rf i - —e ijk a o> 

= -£frk[&& - &#] - i^^ 1 ^ 1 - (104) 

Using !; k = + i£j in Eq. (|104j) and substituting the resulting equation 
into Eq. (JUBJ), we obtain 

- J^t^er - sen - i^^ 7 - 

-^ijM^R - + C^n + CnCT- (105) 

which after a lengthy calculation becomes 

2a 2 M 

m& = - ^Ptf & - ere fi - &?,]• 
-^mme H + -co + zkir + ekf i 

~^a e ^ ajl '- ( 106 ) 

Finally, using for M' the Bondi mass loss equation 

in Eq. (|106jl and replacing M by M = constant , the zero order of the Bondi 
mass we obtain, the equations of motion: 
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9n M 

+f£[&& + - ( 107 ) 

In the last computation we dropped many higher order terms. 
For the Imaginary part, replacing M by M we obtain 



Before we discuss the meaning of Eqs. (|l()7jl and (ll()8j) in the following 
section, we turn to the special case of £j = and a d = <r d . The equations of 
motion then reduce to 

Mo& = %e R + [%c R c R + ^' J ^ir ( 109 ) 

and the Imaginary part 

M a 2 

^ Jk [e R e R } + l^w&ein = o (no) 

which, using Eq. (|l()9j) up to third order, is the identity 

e«*&[M & - 2 ^Cr] = 0. 

As an aside we note that the complex curve in H-space that we wish to 
determine is given by the expression, z a = £ 6 (t), with the complex parameter 
r, while in the above expressions for the determination of the curve we have 
treated £ b to be a complex function of the retarded time Ub, which can take on 
complex values but we have treated as real. In terms of their functional form, 
£ 6 (t) is identical to £ 6 (/Ub), but in terms of the usage in different equations, 
the transition via r = T(ub,C: an d its linearized approximation (done 
in Appendix A) is different and very important. We have, directly, by this 
approximation been throwing out higher order non-linear terms. 
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6 Interpretations, Discussion and Conclusions 



The basic idea we are espousing is that there is, in if-Space, a unique complex 
curve that can be called the intrinsic complex center of mass world-line and 
complex center of charge line. It determines and also is determined from a 
geometric structure on the physical space-time, the null direction field, which 
is given either by the angle field L{ub,C,i or by the tangent field of the 
asymptotically shear- free null congruence l* a . They are geometric quanti- 
ties and exist in the physical space-time independent of the choice of co- 
ordinates or tetrad on 3 + . However the complex curve has no immediate 
physical meaning. There are however a series of observations that do sug- 
gest a physical interpretation. Furthermore there are many clues in earlier 
literatureO EH EH E| EE1 ES| to its significance. 

The basic idea is not to take the parameter space, H-space, as real or 
concrete in any sense - but instead treat it as an observation or holographic 
space. In the physical space-time we are dealing with a complicated physical 
system, a gravitating - charged mass distribution, that creates curvature and 
an electromagnetic field that we can 'see' or observe only from its asymp- 
totic behavior; we can not see individual masses, charges, spins or particle 
trajectories but instead only the gross or large-scale behavior. The idea is 
that the H-space is like a screen that captures certain images of the physical 
space. This idea arises from several directions and, in turn, suggests other 
physical interpretations. 

1. We know how to define energy and momentum from the classical 
treatment of the Bondi mass and momentum. Now from the H-space point 
of view and at a low order of approximation we have recovered the trivial 
relationship between three-momentum, mass and velocity, Eq. ljMl) . i.e., P = 

2. In classical Electromagnetic theory, using the electromagnetic self- 
force and the Lorentz force law, one finds the radiation reaction term f^Cjj ■ 
Here, just by looking at linear order, at the asymptotic Einstein-Maxwell 
fields with no attempt at model building, we obtain exactly the same term in 
the equations of motion without any type of model building. Furthermore, by 
looking at the equations of motion for the special case of £ r = 0, i.e., Eq. (|lfl9|) . 
we see that in addition to the radiation reaction term there are further non- 
linear terms that arises from the Bondi mass loss formula, Eq. (ll03j) . In simple 
numerical examples, neglecting the (a ) 2 term, it appears as if the classical 
run-away solutions of just the radiation reaction term are damped out. 
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3. In our equation for mass loss, Eq. (ll()3j) . there is electric and magnetic 
dipole radiation with the dipole with the second time derivative of the com- 
plex q£, k {uB) = q(£,R + i-e., we see that there is energy loss by both 
electric and magnetic dipole radiation with the identification of q£, k R and 

as the electric and magnetic dipole moments. The numerical factors even 
agree with dipole radiation derived from pure Maxwell theory in Minkowski 
space. 

It appears as if the ^(%) can be used as a 'position vector' and it is 
suggestive that M^(ub) and can be interpreted as the mass and 

charge dipole moments and thus ^r( u b) interpreted as the center of mass 
and charge. 

From knowledge of the Schwarzschild, Kerr, Reissner-Nordstrom and 
charged Kerr metrics, one can try to give meaning to the complex curve 
and in particular to the imaginary part of the curve. It has been well know 
since the 1960's that the Kerr space-time has associated with it a complex 
space and a unique complex world-line such that its imaginary part, = a, 
determines the spin angular-momentum S, via S = Mca. Furthermore for 
the charged Kerr metric there is the same complex world-line with imaginary 
part, |£j| = a from which both the spin S and magnetic dipole moment, /x, 
are obtained from S = Mca and \i = qa. In none of these cases is there a 
non-vanishing real part of the world-line and hence the real center of mass 
and the electric dipole vanished and, hence, are at the origin. 

4. These observations suggest that in the present work the imaginary 
part of the complex curve, i.e., £ f , should be used to define both the spin 
and the magnetic dipole moment via 

~S = Met I 

5. From a stndv [T8llT2lllT] of the Maxwell equations in Minkowski space, 
where one needs far less approximations and sometimes no approximations 
at all, again we have that ~Ji = q £ /. 

6. Though at the present time we do not know the meaning of many of 
the terms in both the expression for the 3-momentum, (|9"4"J) and the equations 
of motion, (|107[) . it appears possibly that they come from the fact that we are 
dealing with the matter from entire interior of the space-time and information 
from earlier times is being fed in via the higher time derivatives. Nevertheless 
some of the terms are familiar: In (|94j) there are three terms 
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p k = m& - - y c m ^^r + 

that are known from other considerations: there is the standard Mv term, 
the second term comes from radiation reaction [its time derivative is the well- 
known radiation reaction force.] The last term is the spin-velocity coupling 
term found in the Mathisson-Papapetrou equations[2Il[|[2TJ if again we identify 
the spin with Mc^j. 

It thus appears natural, in our case where our two world-lines coincide, 
to identify £j as the imaginary center of mass and charge and thus define the 
full complex ^(ub) as the complex center of mass and center of charge. This 
point of view allows a unification of many examples of electric and magnetic 
phenomena. 

Returning to the charged Kerr case, it was noted that in the ratio S/fM — 
Mc/q, the imaginary displacement, 'a', drops out, and one discovers that for 
the charged Kerr solution the gyromagnetic ratio has the Dirac value, i.e., 
g = 2. This result came from the fact that the two complex world-lines 
coincided. From this observation we see that, for all asymptotically flat 
Einstein-Maxwell fields where the two complex H-space world-lines coincide, 
again the gyromagnetic ratio will have the Dirac value, g — 2. 

Though there is still much to understand and do, we believe that there 
is something essentially correct in this screen or holographic point of view 
towards these H-space curves. A variety of questions, however, do arise. 

a. Is this discovery of the complex curve on H-space and the direction 
fields just a booking keeping device or is there something deeper. We do not 
know but we can see a beautiful symmetry between orbital and spin angular 
momentum and between electric and magnetic dipole moments. They are 
naturally unified into the complex vector £ a . This structure takes on even 
more beauty in the case of the twisting type II metrics (in some sense funda- 
mental solutions) where the asymptotically shear-free congruence is totally 
shear-free. 

b. We do not yet know the connection between our definition of spin 
and the more usual one obtained from the symmetry properties of the BMS 
group. 

c. We need to analyze the meaning of many of the terms in the equa- 
tions of motion. It appears clear that they will involve couplings between 
all the different types of moments; mass-dipole/orbital angular momentum, 
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spin, electric and magnetic dipoles and higher moments. Will they have any 
obvious meaning? 

d. There is the issue to be explored of whether it is possible in principle 
[almost certainly not in practice] to observe the complex world-line. The 
answer very likely is yes; it probably can be done by a spherical harmonic 
analysis of the angle field, L(ub, (, (), then looking at the I = 1 component. 

e. The potential damping of the classical run-away solution from the 
radiation reaction term must be further studied. 

f. What relationship does our unique choice of the angle field L(ub,(X) 
have to the same geometric object in the algebraically special metrics where it 
is used to define a CR structure on 3 + ? It would appear that a mathematical 
statement of our results could be the following: (i). For all asymptotically 
flat vacuum space-times there is a unique CR structure on ii). For 
all asymptotically flat Einstein-Maxwell space-times there are two unique 
CR structures on 3 + , (iii). In the degenerate case, when the two structures 
coincide, one find that the gyromagnetic ratio is that of Dirac, namely g — 2. 
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8 Appendix A 

In order to obtain r = T(ub, C> C)> the inversion of 

u B = X(t, C, C) = C(r)L(C, C) + X l > 2 {r, c, a 

we first note that r can be replaced by any function of r, i.e., the Eq. (|l9|) is 
invariant under r = F(t). See Eq.fJHHI). Thus by taking 

e°(r) = V2r 

we have 

^ = r + r(r)/ l (C,C)+^>2(r,C,C) 

or 

r = u B - C(rMC,C) + X^tXX)- 
This can be solved by iteration to any order: 

To = u B 

ri = u B -e(u B )k(cx)+xi> 2 (ux,c) 

T n = U B ~ C(T n -l)k((,() +Xj> 2 (r n _i,C,C)- 

In the text we have truncated this procedure and used the approximate 
inversion as 

T = U B - C(u B )k(CX)- 

9 Appendix B 

Normalization and conventions: Our goal here is to express the function \1/ in 
terms of the Bondi mass/momenta and also show how to extract them from 
the 

We start, in a given frame in Minkowski space, with C a a unit radial 
space-like vector and T a a unit time-like vector given by 

C a = (0, cos^sintf, sin0sin6>, cos#) = —C a , (111) 
T a = (1,0,0,0). (112) 
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with the null vectors 

l a = —i=(l, — cos(f)sm9, — sm(f)sm9, — cos9) — —i=(T a + C a )(113) 
v 2 v 2 

n a = —= (l,cos<f)sm6,sm(f)sm6,cos6) = —=(T a — C a ) (H4) 
v2 v2 

and 

c a = ^ a — n a — —V2(0, cos 4> shi 0, sin sin 8, cos 0) = V2C a 

We let a = (0, 1, 2, 3) = (0, i) and / a = (l , k) and q = V2Q, etc. 
We then define the / = and I = 1 parts of the mass aspect \& by 



V = ip* 2 ° + 2LQW +L 2 a" + tfa + aa (115) 
* = + (116) 

Therefore 

Ma = X/a-X i Q/ a + ... (117) 
= [ X T a - QC a + X C a - V^x'CiTa] 



From the integral identities, with dS the area element on the unit metric 
sphere, 

dS = 4tt 



/ 



CAS = 



Aix 

CiC jdS = „ Sn 
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we have that 



J m a ds = -±= J dS( X T a - V2 X l c t c a ) 



4ir 



V2 



V2 



V2 



(118) 
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with Xi = -X 1 - 



On the other hand, from reference and a change of notation [see 
footnote 1], we have that the 4-momentum for an asymptotically flat space- 
time is defined by 



c 3 



P a = (Mc,-P l ) = -— m a dS, (119) 

07TCx J 

C £\C 

* = -2V2-M + — P% + %> 2 (120) 



then from Ens. (fTT9l) . (fT2nl and fTT7jl . 



(P ,P,) = (M C) P,) = -^^ (x>**Y 
We find that: 

X = -2V2-M (122) 
X, = 6^P\ (123) 

10 Appendix C 

Relations between the following quantities have been used through this work: 



" 2 u 'i+cr i+cr i+cr' ( j 



2 i + cc i + CC l + CC 



2 i + CC l + CC i + CC 
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t a = ^(1,0,0,0), (125) 

2 i+_CC i + CC i + CC _ 

c * = i-_ n . =v ^( 0) i±i,-<-i^i,zi±^). 

l + CC i + CC i + CC 

Letting a,b, etc., take the values (0, 1,2,3) = (0, i) we have the products 
( - - \ 



(rriiCj - rrijCi) = -iV2e ijk m k 

1/ \ 1 

-{rrncj - m jCi ) + - 

-i—=eijkin k + (I = 2)Terms 



rrii cj = -{rriiCj - rrij d) + -{rriiCj + m.,- Cj j 



V2 

In a recent article [TH] we have developed a notation to describe symmetric 
and trace-free tensor products of the three Euclidean vectors (c^m^mj) and 
their products; 

Several examples that have been used in the text are: 

C - - -r (0) 
m = Y$ 

m i™J = Y (2)ij 
CjTflj + mjCj = — Y^2)ij 
•\rr —9fi- — 

where the upper index indicates the spin weight and the lower index in paren- 
thesis gives the I value. The indices i and i,j, etc., are Euclidean tensor 
indices. 

Some examples of product decompositions (Clebsch-Gordon expansion) 
that have been used in the text are: 

ylyO _ * - V l , V 1 
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We have found that this type of notation and analysis is extremely useful. 
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